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Abstract 

Neutrino masses may be generated by the VEV of an SU{2)i Higgs triplet. We 
assume that the doubly charged component of such a triplet has a mass in the range of 
several 100 GeV, such that it is accessible at LHC. Its decay into like-sign leptons pro- 
vides a clean experimental signature, which allows for a direct test of the neutrino mass 
matrix. By exploring the branching ratios of this decay into leptons of various flavours, 
we show that within this model the type of the neutrino mass spectrum (normal, in- 
verted or quasi-degenerate) might actually be resolved at the LHC. Furthermore, we 
show that within the Higgs triplet model for neutrino mass the decays of the doubly 
charged scalar into like-sign lepton pairs at the LHC provide a possibility to determine 
the Majorana CP phases of the lepton mixing matrix. 
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1 Introduction 



Recent developments in neutrino physics demand for an extension of the Standard Model in 
order to give mass to the neutrinos. A popular way to achieve this goal is to introduce right- 
handed singlet neutrinos. An alternative, equally valid and rather economical possibility 
is to extend the scalar sector of the Standard Model. In addition to the Higgs doublet, 
scalar representations consistent with the SU{2)l x U{1)y gauge group and the possible 
fermionic bilinears are pQ a triplet, a singlet with charge +1, or a singlet with charge +2, 
see Refs. [21 [3l H] for corresponding models. In this work we focus on the first mentioned 
possibility, namely a scalar SU{2)l triplet. Such a triplet arises naturally in many extensions 
of the Standard Model, for example in left-right symmetric models [5], or in Little Higgs 
theories [U [7j. When the neutral component of the triplet acquires a vacuum expectation 
value (VEV), vt, a Majorana mass term for neutrinos is generated at tree level, proportional 
to Vt- In order to obtain small neutrino masses this VEV and/or the corresponding Yukawa 
couplings have to be very small. If a very high energy scale M 3> f = 246 GeV is associated 
to the triplet, one obtains the well-known seesaw (type-II) relation vt ~ v'^/M as explanation 
for the smallness of neutrino masses [HI El [IH] ■ 

Here we consider a different scenario, assuming that the triplet states have masses not 
too far from the electroweak scale. The Higgs potential of the Standard Model Higgs doublet 
(j) and the triplet A contains a term /i0A0, which breaks lepton number explicitly. Assuming 
that all other mass parameters in the potential are of the electroweak scale v, the minimi- 
sation of the potential leads to the relation for the triplet VEV vt ~ /x, see e.g. [H]. The 
hierarchy fi <^ v may find an explanation for example through extra dimensions [12]. A 
Higgs triplet slightly below the TeV scale is the generic situation in Little Higgs theories [6], 
see Ref. [7] for a discussion of neutrino masses in this framework. Other examples for models 
with TeV scale triplets responsible for neutrino masses can be found, e.g., in Refs. [131 [H [IS] • 
Our phenomenological analysis does not rely on a specific model realisation, apart from the 
assumption that neutrino masses arise from a triplet with masses in the TeV range. 

The hypothesis of such a Higgs triplet can be tested at collider experiments. In particular, 
if kinematically accessible, the doubly charged component of the triplet will be produced 
in high energy collisions, and its decay into two equally charged leptons provides a rather 
spectacular signature, basically free of any Standard Model background. This process has 
been studied extensively in the literature (see Refs. [T6l[I7l[IHl[l9l[2Ql[2Il[22l[23ll21] for 
an incomplete list), and has been used to look for doubly charged scalars at LEP [23] and 
Tevatron [26]. These searches resulted in lower bounds for the mass of the order Mh++ ^ 
130 GeV. Therefore, we will consider in the following masses in the range 130 GeV < Mh++ < 
1 TeV, above the present bound but still in reach for LHC. 
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If the Higgs triplet is responsible for the neutrino mass the decay rate for H^^ — i'^i'l 



is proportional to the modulus of the corresponding element of the neutrino mass matrix 
|Mabp. This opens a phenomeno logically very interesting link between neutrino and collider 



physicqj, and by the observation of like-sign lepton events at LHC a direct test of the neutrino 
mass matrix becomes possible. In this work we assume that a doubly charged Higgs is 
indeed discovered at LHC, and we use the information from the decays H^^ —>■ i^K ^'^ 
learn something about neutrinos, under the hypothesis that the neutrino mass matrix is 
dominantly generated by the triplet VEV. 

Current neutrino data leave some ambiguities for the neutrino mass spectrum. The 
neutrino mass states can be ordered normally or inverted, and the masses can be hierarchical 
or quasi-degenerate. We will show that under the above assumptions actually LHC might 
play a decisive role in distinguishing these possibilities. Furthermore, we show that it might 
be possible to determine the Majorana phases P, |33] in the lepton mixing matrix, which 
in general is a very difficult task. Implications of the different possibilities of the neutrino 
mass spectrum for the decay of a doubly charged scalar in the Higgs triplet model have been 
considered previously in Ref. [M], see also p2]. Building upon the results obtained there, 
we perform a full parameter scan including all complex phases, which — as we will see — play 
a crucial role for the relevant observables. 

The outline of the paper is as follows. In Sec. [2] we present the general framework, where 
in Sec. 12.11 we review how the neutrino mass matrix arises in the Higgs triplet model, and in 
Sec. 12.21 we discuss the signature of the model at LHC. Sec. O contains the main results of 
our work. After describing our analysis in Sec. 13.11 we discuss in Sec. 13.21 how the branching 
ratios of the doubly charged scalar depend on the parameters of the neutrino mass matrix. 
In Sec. I3.3l we investigate the possibility to determine the type of the neutrino mass spectrum 
from like-sign lepton events at LHC, whereas in Sec. 13.41 we show that within this framework 
indeed Majorana phases can be determined. Concluding remarks follow in Sec. HI 

2 Framework 

2.1 Neutrino masses from a Higgs triplet 

If an SU (2) l Higgs triplet with hypercharge F = 2 is present in the theory the following 
renormalisable term appears in the Yukawa sector of the Lagrangian: 



where the indices a,b = e, jjL,T label flavours. La are the lepton doublets, C is the charge 
conjugation matrix, T2 is the Pauli matrix, A denotes the scalar triplet, and fab is a symmetric 

^Such a link exists also in other classes of models, see for example j271 [28l [29l [30l [3T1 [32] , However, in 
most cases the connection between collider signals and the neutrino mass matrix is much less direct as in 
the Higgs triplet model. 




(1) 
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complex Yukawa matrix. Without loss of generality we work in the mass basis of the charged 
leptons. The components of the triplet are given by: 

A - ( ""'^ \ (2) 

The VEV of the neutral component {H'^) = Vt/\/^ induces a Majorana mass term for the 
neutrinos: ^ 

-vl,C-^Mal>VLb + ^.C. with Mab = V2VTfal>- (3) 

We assume in the following that this is the sole source for neutrino masses (or at least the 
dominant contribution). As usual the neutrino mass matrix Mab is diagonalised by: 

M = ?7diag(mi, m2,ms)U'^ . (4) 

For the PMNS matrix U we adopt the parametrisation 



U = V diag(e* ^ ,e' ^ , e*^) with 



V = 



~C23Si2 — Si3S23Ci2e*'^ C23C12 — Si3S23Si2e*'^ S23C13 
\ S23S12 - Si3C23Ci2e*'^ -S23C12 - Si3C23Si2e*'' C23C13 / 



(5) 



where 6 is the so-called Dirac CP violating phase which is in principle measurable in neutrino 
oscillation experiments, and are the Majorana phases [HI [33]. Note that only relative 
phases physical, and therefore there are only two independent Majorana 

phases. Neutrino oscillation data determine the so-called solar and atmospheric oscillation 
parameters [35] : 

sin^ = 0.32 ± 0.023 , Amj^ = (7.6 ± 0.20) x 10^^ eV^ , 

(6) 

sin^^ss = 0.50 ±0.063, |Am^J = (2.4 ±0.15) x 10~^eV^ 

where we give la errors and Amfj = mf — m'j. For the mixing angle ^13 there is only an 
upper bound, 

sin^^i3<0.05 at 3a, (7) 

whereas nothing is known about the phases 6,aij. The ordering of the mass states is de- 
termined by the sign of Amg^: for normal hierarchy (NH) Am^^ > 0, whereas for inverted 
hierarchy (IH) we have Am|]^ < 0. We denote the lightest neutrino mass by mo, hence, 

'»» - f am ■ W 

[ m3 (IH) 

If mo > a/I Amgj^l ~ 0.05 eV the neutrino mass spectrum is quasi-degenerate (QD). The 
most stringent bound on the absolute scale of the neutrino mass comes from cosmology, 
which is sensitive to the sum of the three masses. In a recent analysis ^36j the upper bound 
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Ylii'^i < 0-5 eV at 95% CL has been obtained, which translates into mo < 0.16eV. Since 
this corresponds to the QD regime the bound is the same for NH and IH. Taking into account 
Eq. ([3]), the constraint from cosmology applies directly to the product of triplet VEV and 
Yukawas: 

yT/a;>< 10-^° GeV. (9) 
2.2 Doubly charged scalars at the LHC 

At the LHC the process 

pp H++H— rr (10) 

provides a very spectacular signature, namely two like-sign lepton pairs with the same invari- 
ant mass and no missing transverse momentum, which has essentially no Standard Model 
background. The pair production of the doubly charged scalar occurs by the Drell-Yan 
process qq — ^ 7*,Z* H H^^ ^ with a sub-dominant contribution also from two-photon 
fusion 77 H H^^ . The cross section is not suppressed by any small quantity (such 
as the Yukawas or the triplet VEV) and depends only on the mass Mh++, see e.g. [IHl [21] • 
QCD corrections at next-to-leading order have been calculated The cross section for 

H H^^ pair production at the LHC ranges from 100 fb for a Higgs mass Mh++ = 200 GeV 
to 0.1 fb for Mh++ = 900 GeV [2lj. Hence, if the doubly charged scalar is not too heavy a 
considerable number of them will be produced at LHC assuming an integrated luminosity of 
order 100 fb'^ 

The rate for the decay if — > is given by 

T{H++ - W) = ■ , \fa,?MH.. , (11) 

47r(l + dab) 

with Sab = 1(0) foT a = b{a ^ b). Hence, the rate is proportional to the corresponding 
element of the neutrino mass matrix |Mabp. This observation is the basis of our analysis. 
Using Eqs. and fill I) the branching ratio can be expressed as 



and from Eq. (jll) and the unitarity of U follows 

\^IM |2 _ V- ^2 _ / Smg + Ami + A^li (NH) 
2^171/4,1 -2^ m, - <i g^,^ ^^2^^ 2|Am2j (IH) ' ^^^^ 



cd i=l 



In addition to the lepton channel the doubly charged Higgs can in principle decay also 
into the following two-body final states including singly charged Higges and/or the W: 

H++ , if ++ ^ , ii++ ^ W^W-^ . (14) 



^Let us note that — depending on the mass sphtting between the double and single charged components 
of the triplet — also the channel q'q H^^H^ may significantly contribute to the production of doubly 
charged scalars, see e.g. [THI12]. 
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The first two decay modes depend on the mass sphtting within the triplet. We assume in 
the following that they are kinematically suppressed. The rate for the WW mode is given 
by 



r(i/++ ^ W+W+) ^ ^-^iP^ , (15) 

where v = 246 GeV is the VEV of the Standard Model Higgs doublet, and we have used 
Mh++ 3> see e.g., Ref. [21] for full expressions and a discussion of possibilities to 

observe this process at LHC. Hence, the branching ratio between t^t^ and W^W^ decays 
is controlled by the relative magnitude of the triplet Yukawas fab and the VEV vt- The 
requirement r(iJ++ W^W^) < T{H^^ — > ^q^^), together with the constraint from 
Eq. ([9]) implies: 

^x<i,-o(m^\'\ (16) 



Mh++ J 

The triplet VEV contributes to the p parameter at tree level as |2] p ~ 1 — 2{vj'/vY. 
The constraint from electroweak precision data p = 1.0002l°:°°g4 

at 2a [37] translates into 

vt/v < 0.02, which is savely satisfied by requiring Eq. ffTBj) . 

In this model contributions to lepton flavour violating processes, — 2, and in principle 
also to the electron electric dipole moment are expected, see e.g. [HH [HHl [39] and references 
therein. Following Refs. [HI [39], the most stringent constraint on the Yukawa couplings fab 
comes from p, —>■ eee, a process which occurs at tree level via Eq. ([T]). The branching ratio 
for this decay is given by [39] : 

BR(, ^ eee) . « 20 [-^) " IfM^ . (17) 

Hence, the experimental bound BR(/i — > eee) < 10~^^ [37J constrains the combination 
\feefeti\ ^ 2 X 10"'^(Mi^++/100 GeV)^. Assuming that all fab have roughly the same order of 
magnitude we obtain an estimate for the interesting range of the Yukawa couplings: 



where the lower bound emerges from Eq. ( 11611 assuming that the bound ([9]) is saturated. We 
see that several orders of magnitude are available for the Yukawa couplings. For fab close 
to the lower bound of Eq. f[T^ the decay H~^^ W^W^ will become observable at LHC, 
whereas close to the upper bound a signal in future searches for lepton flavour violation is 
expected, where the details depend on the structure of the neutrino mass matrix [Sll [32] • 
The interval for the Yukawas from Eq. ( fTSi) implies a triplet VEV roughly in the keV to MeV 
range. 

The basic assumption in our analysis is that a sufficient number of like-sign leptons is 
observed. If some of the decay modes of Eq. (IT^ are present the number of dilepton events 
will be reduced according to the branching. If enough events from both types of decay 
(leptonic and non-leptonic) were observed in principle an order of magnitude estimate for 



6 



the Yukawa couplings fab and the triplet VEV vt might be possible [HI 122] • Here we do 
not consider this case and use only dilepton events, and therefore, we do not obtain any 
information on the overall scale of the fab in addition to Eq. (ITSll . 

3 Numerical analysis and results 
3.1 Description of the analysis 

As mentioned above, we focus in our analysis on the process fllOl] . which provides the clean 
signal of four leptons, where the like-sign lepton pairs have the same invariant mass, namely 
the mass of the doubly charged Higgs. Given the fact that the branching H~^^ ^t^t 
proportional to the neutrino mass matrix, one expects all possible flavour combinations of 
the four leptons to occur, including lepton flavour violating ones. In Ref. [21] simple cuts 
have been defined for final states consisting of electrons and muons, eliminating essentially 
any Standard Model background. 

In general tau reconstruction is experimentally more difficult because of the missing trans- 
verse energy from neutrinos. However, in the case of interest enough kinematic constraints 
should be available to identify also events involving taus. It turns out that the inclusion 
of such events significantly increases the sensitivity for neutrino parameters. Therefore, fol- 
lowing Ref. [^, we assume that events where one of the four leptons is a tau can also be 
reconstructedo This should be possible efficiently, despite the complications involving the 
tau reconstruction, since the invariant mass is known from decays without tau, which can 
be used as kinematic constraint for events of the type C^^t^ for i = e or fi. Further- 
more, one can adopt the assumption that the neutrinos carrying away the missing energy 
are aligned with the tau. 

In principle it is difficult to distinguish a primary electron or muon from the ones orig- 
inating from leptonic tau decays. Since here we are interested in investigating the flavour 
structure of the decays, leptonically decaying taus might be a "background" for the Higgs 
decays into electrons and muons, and vice versa. However, due to the energy carried away 
by the two neutrinos from the leptonic tau decay, a cut on the invariant mass of the like-sign 
leptons should eliminate such a confusion very efficiently. It is beyond the scope of this 
work to perform a detailed simulation and event reconstruction study. The above arguments 
suggest that our assumptions are suitable to estimate the sensitivity of the Higgs decays for 
neutrino parameters by the procedure outlined in the following. 

We define as our five observables the number of like-sign lepton pairs with the fiavour 
combinations 

x = (ee), (e/i), (/i/i), (er), (/ir) . (19) 
Note that these five branchings contain the full information, since BRt-,-, which we do not use 

•^To be conservative we do not include events with more than one tau, since aheady the inclusion of events 
with one tau provides enough information for our purposes. 
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explicitly, is fixed by BR,-,- = 1 — Ylx^^^- Taking into account the number of occurrences 
of the combinations ffTI?|) in four leptons where at most one tau is allowed, the number of 
events in each channel is obtained as: 

Nab = 2N2H e BRab ^ BR^ for (ab) = (ee), (e/i), {nfx) , 

(20) 

Nab = 2N2H e BR,,(BRee + BR^^ + BR^^) for (ab) = (er), (/xr) , 

where N2H is the total number of doubly charged scalar pairs decaying into four leptons, 
and e is the detection efficiency for the four lepton events. For simplicity we assume here 
a flavour independent efficiency. The branching ratios are given in Eq. ( fT2l) . To illustrate 
the sensitivity to neutrino parameters we will use eN2H = 10^ or eN2H = 10^ events. For an 
integrated luminosity of 100 fb~^ at LHC these event numbers will be roughly obtained for 
Mh++ ^ 350 GeV and Mh++ ^ 600 GeV, respectively [21]. 

To carry out the analysis we define a function from the observables in Eq. (|20|) . 
For given eN2H they depend only on neutrino parameters. We consider five continuous 
parameters: the lightest neutrino mass mo, S13, the Dirac phase 6, and the two Majorana 
phases ai2 = ai — 0:2 and 0:32 = 0:3 — 02, plus the discrete parameter h = NH or IH describing 
the mass ordering. The remaining neutrino parameters, the two mass-squared differences 
and the mixing angles S12 and S23, are fixed to their experimental best fit values given in 
Eq. (E]). The is constructed as: 

x'^imo, si3, 5, ai2, 032, ^) = K S^^ Vy + — with 

WJ (21) 

K = N^'^imo, si3, <5, 012, a32, h) - iVf p 

where x and y run over the five combinations given in Eq. ffTIJl) . For the "data" iV^^^ we use 
the prediction for A^^; at some assumed "true values" of the parameters, 
(mo, Si3, (5, 012, 032, ^)*™°- Then the statistical analysis tells us the ability to reconstruct 
these true values from the data. For the covariance matrix S we assume the following form: 

S^y = N^d^y + aL^iVP-^iVP-^ + . (22) 

It includes statistical errors, a fully correlated normalisation error o"nomi, and the uncertainty 
introduced from the errors on the oscillation parameters 5'°'^'^. The normalisation error o"norm 
arises from the uncertainty on the luminosity and the efficiency. Moreover, the possibility 
that the non-leptonic decays of H^^ of Eq. ([T3j) might occur at a sub-leading level and 
are not observed introduces an uncertainty in the number of leptonic decays. We adopt a 
value of (Tnorm = 20%. We have checked that even an analysis with free normalization (i.e., 
Cnorm oo) Icads to Very similar results. This means that the information is fully captured 
by the ratios of branchings^ 



''This is true as long as all branchings from Eq. are used; if the events containing taus are omitted 
our results depend to some degree on the value adopted for tTnorm- 
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Via the covariance matrix 5°^"^ we account for the fact that the parameters Am^i, |Am|^|, 
Si2 and S23 have a finite uncertainty. We include the errors from Eq. and take into account 
the correlations which they introduce between the observables N^. The last term in Eq. (|2T1) 
takes into account the constraint on S13 from present data according to Eq. ([7]). Let us note 
that within the time scale of a few years the errors on oscillation parameters are likely to 
decrease. In particular, also the bound on S13 will be strengthened or eventually a finite 
value could be discovered by upcoming reactor or accelerator experiments, see for example 
Ref. [lO]. To be conservative we include only present information, although at the time of the 
analysis better constraints might be available. We have checked that the precise value of S13 
within the current limits as well as its uncertainty have a very small impact on our results, 
and a better determination may lead at most to a marginal improvement of the sensitivities. 

3.2 Branching ratios 

In Fig. [1] we show the branching ratios for NH and IH as a function of the lightest neutrino 
mass mo- For fixed mo, the interval for the branching emerges due to the dependence on 
the phases ai2,a32,S, and also the uncertainty on solar and atmospheric oscillation param- 
eters contributes to the interval. In the plots one can identify the regions of hierarchical 
neutrino masses, mo < 10~^ eV, and QD masses, mo > 0.1 eV, where NH and IH become 
indistinguishable. In the limiting cases mo = and mo 00 the analytic expressions for 
the branchings are rather simple. For NH and mo = one finds to leading order in the small 



quantities r = Amg^/jAmg^l ^ 0.03 and s^g < 0.05 (at 3cr): 

Bj^NH,™o=o ^ 42^ + 244v^cos(a32-25), (23) 

gj^NH,™o=0 ^ 2[s?244r + s234 + 2si2Ci2S23C23Sl3V^COs(a32-5)] , (24) 
Bj^NH,r„o=0 ^ 4 + 2444v^COSa32 + 44r 

-4s23C23Sl2Cl2Si3A/r cos(a32 - S) , (25) 

Bj^NH,mo=0 ^ 2[s?2C?24^ + 44-2Si2Ci2S23C23Si3V^COs(«32-^)] , (26) 

Bj^NH,™o=o ^ 2sl,cl, (1 - 2cLv^cos a32 + 4r) . (27) 
For IH and mo = 0, S13 = we have 

BR« = l(l-sin2 2e,2W^) , (28) 

Bj^iH,™o=o ^ sin' 26,2 sin' ^, (29) 

BRf-o^o = ^ (1 - sin^ 26,2 sin' ^) , (30) 

BR^»'-"=° = 4sin2 2^i2sin2^, (31) 

BRf^'-o=o ^ 1 3i^2 2^^3 _ gi^2 2^^^ ^.^2 ^ j ^ ^32) 
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Figure 1: Branching ratios BR(_ff (-a^b) as function of the hghtest neutrino mass toq for NH (hght-red) 
and IH (dark-blue). The thick sohd hnes are for s\z — 0, and the thick dashed hnes for S13 — 0.1, where 
the dependence on phases as weh as the uncertainty of solar and atmospheric oscillation parameters at 2cr 
are included. The thin solid lines show the branchings for oscillation parameters fixed at the best fit points 
Eq. ([U, Si3 = 0, a-ii = tt, and ai2 = 0, 7r/4, n/2, 3n/A, n. 



and in tlie limit mo — ^ 00 with S13 = the branchings become 
BR?P = sin^ 2eu sin^ ^) = ^ bR^»'-o=° , 



brQ° 



C23 sin^ 26*12 sin 



2 



1 — - sin^ 26*23 (1 ■3^2 cos 0:31 — c^2 COS 0:32) — C23 sin^ 26*12 sin^ 
S23 sm 26^12 sm — - - BR^^' , 



(33) 
(34) 
(35) 
(36) 



^ sin^ 26*23 f 1 - ■S12 cos 0:31 - c^2 cos 032 - ^ sin^ 26*12 sin^ ^ ) . (37) 



Note that for a vanishing hghtest neutrino mass, mo = 0, there is only one physical Majorana 
phase, 032 for NH, and 0:12 for IH, as clear from Eqs. (jl]) and ([5]). 
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In the following we will explore the parameter dependencies of these branchings to obtain 
information on the neutrino mass spectrum and on Majorana phases. The rather wide ranges 
for the branchings in the cases of IH and QD spectrum suggest a strong dependence on the 
phases, and as we will see in Sec. 13.41 these are the cases where Majorana phases can be 
measured very efficiently. The determination of the mass spectrum is somewhat more subtle. 

A clear signature for the NH with small mo is provided by BRefEI. Eq. ( l23l) shows that 
for NH and mo = 0, BRge is suppressed by r and/or sf^, and there is the upper bound 
BRee < 5.3 X 10^^ for the largest value of allowed at 2a and s^g = 0.01, in agreement 
with Fig. [H In contrast, for IH with mo < 0.01 eV and for QD spectrum, Eqs. ( l28l) and ( l33l) 
give the lower bounds BR^^ > (1 - sin2 2^i2)/2 ^ 0.03 and BRee > (1 - sin2 2^i2)/3 ^ 0.02, 
respectively. Therefore, the characteristic signature of normal hierarchical spectrum is the 
suppression of Higgs decays into two electrons. 

From a first glance at Fig. [T] one could expect that it might be difficult to distinguish IH 
and QD spectra, since there is always overlap between the allowed regions in the branchings. 
Indeed, if only branchings involving electrons and muons (BRee, BRe^, BR^^) are considered 
there is some degeneracy between IH and QD, especially if S13 is allowed to be close to the 
present bound. However, as we will show, due to the complementary dependence on the 
phases of all the BKab including also taus, the degeneracy is broken and these two cases can 
be disentangled. Consider, for example, BR^^ and BR^,-: in the case of IH with mo = 
they behave very similar as a function of ai2, see Eqs. (!30|) and (132|) . whereas for QD they 
show opposite dependence, compare Eqs. flH3|) and flHTj) . and phases which give BR|^° = 
maximise 

Note that for S13 = and S23 = 0.5, BRe^ and BRer are identical. Nevertheless there 
is important complementariness between them. First, the uncertainty on S23, see Eq. dH]), 
affects each of them significantly, and it reduces the final sensitivity if only BRe^ is used 
in the analysis. But since BRe^ and BRe,- are related by the transformation S23 —>■ C23, 
C23 —S23 this uncertainty is cancelled if both of them are included in the fit. Second, 
it can be shown that the leading order term in S13 is the same for BRe^ and BRer, apart 
from an opposite sign. Therefore, also the impact of S13 is strongly reduced if information 
from both of them is taken into account. One can observe from Fig. [1] that for small mo 
and NH, BRe^ and BRe,- show a significant dependence on S13, while in the other cases the 
dependence is mild. The reason is a leading term linear in ^/rsl3 in Eqs. ( 12^ and ( 126|1 . 
whereas in all other cases S13 appears either in sub-leading terms or at least at second order. 

3.3 Determination of the neutrino mass spectrum 

Let us now quantify the ability to determine the neutrino mass spectrum by performing a 
analysis as described in Sec. 13.11 In Fig. [2] we show the by assuming that "data" are 
generated by a hierarchical spectrum with normal ordering (left), a hierarchical spectrum 

^Note that the behaviour of BRee is the same as the effective neutrino mass probed in neutrino-less double 
beta-decay, which is also proportional to |Mee|, see for example Ref. |41) . 
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10'^ 10"^ 10"^ 10"^ 10° 10"^ 10"^ 10"^ 10° 10"^ 10"^ 10"^ io' 
m, (eV) m, (eV) m. (eV) 



Figure 2: Xmin '^o assuming a true hierarchical spectrum with NH (left) and IH (middle), and a true 
QD spectrum (right). The is shown for eN2H = 100 (dashed) and 1000 (solid) events, and (Tnorm — 20%. 
In the fit we assume either NH (light-red) or IH (dark-blue), and we minimise with respect to S13 and the 
phases. We adopt the following true parameter values. Left: mg = 0, NH, 032 = tt; middle: mo — 0, IH, 
ai2 = 0; right: too = 0.15 eV, ai2 = O.Itt, 0:32 = I.Gtt; and always S13 = 0. 

with inverted ordering (middle), or a QD spectrum (right). These data are fitted with both 
possibihties for the ordering (NH, light-red curves, and IH, dark-blue curves) and a value for 
mo shown on the horizontal axis. We minimise the with respect to the other parameters, 
taking into account the current bound on S13. The results are shown for a total number of 
doubly charged scalars decaying into like-sign leptons of eN2H = 10^ (solid) and 10^ (dashed). 

First we discuss the sensitivity to hierarchical spectra with a very small lightest neutrino 
mass mo . The left panel of Fig. [2] shows that a NH with small mo can be identified with very 
high significance. An inverted hierarchical spectrum as well as a QD spectrum have Ax^ ^ 60 
already for 100 events. An upper bound on the lightest neutrino mass of mo < 0.01 eV at 
3a can be established by LHC data. As discussed in the previous section this information 
comes mainly from the suppression of the decay into two electrons, which occurs only for 
normal hierarchical spectrum. An inverted hierarchical spectrum (middle panel) can be 
distinguished from a QD one at around 3cr with 100 events, where the increases roughly 
linearly with the number of events. The ability to exclude a QD spectrum in case of a true 
IH depends on the true value of the Majorana phase au- The example chosen in Fig. [2l 
Q/*^"^ = 0, corresponds to the worst case; for all other values of the for QD is bigger. 

Fig. [3] shows the ability to identify a hierarchical spectrum as a function of the true value 
for the Majorana phase, where for mo = there is only one physical phase. The shaded 
regions show that for 1000 events the true spectrum can be identified at 5a significance, and 
an upper bound on the lightest neutrino mass mo < 8 x 10^'^ eV for NH and mo < 4 x 10~^ eV 
for IH is obtained, independent of the true phase. For the black contours in Fig. [3] we do not 
use the information from decays into taus, i.e., we use only the lepton pairs (ee), (e/i), (/i/i). 
This analysis illustrates the importance of the tau events. For example, if tau events are 
not used an IH with mo = cannot be distinguished from a QD spectrum for a^™'' ~ vr. 
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Figure 3: Determination of hierarchical neutrino mass spectra, mg™ — 0, assuming 1000 Higgs pair 
decays. The upper (lower) panels are for a true NH (IH), and for the left (right) panels the fit is performed 
assuming a NH (IH). As a function of the true value of the Majorana phases we show contours — 4, 9, 16, 25 
(from dark to light), minimising with respect to all parameters except from mo. Coloured regions correspond 
to our standard analysis, whereas for the black contours we do not use decays into tau leptons. 



Also the sensitivity to a NH is significantly reduced, which becomes even more severe if less 
events were available. 

Now we move to the discussion of a true QD spectrum. As shown in the right panel of 
Fig. [2] also a QD spectrum can be identified quite well, and a lower bound on the lightest 
neutrino mass of mo > 2 (6) x 10^^ eV at 3cr can be obtained for 100 (1000) events. Note 
that for the example shown in Fig. [21 100 events give a A^^ ^ 12.4 for the IH with tuq = 0, 
which corresponds roughly to an exclusion at 3.5a. The potential to exclude a hierarchical 
inverted spectrum depends on the true values of the Majorana phases, and the true values 
of ai2 and 0:32 adopted in Fig. [2] correspond to the worst sensitivity. In Fig. H] we show 
contours of Ax^ for IH with mo = assuming a true QD spectrum, in the plane of the true 
Majorana phases. For 1000 events we find some islands in the plane of 0:12 and 0:32 where 
the reaches values as low as 30 (compare Fig. [2]), however in most parts of the parameter 
space the exclusion is at more than 7a. For 100 events typically a significance better than 
4cr is reached, but there are some notable regions (— 7r/2 < 0:12 < vr/2 and 032 ~ 7r/2, 37r/2) 
with values between 16 and 9. 

Let us add that for the exclusion of an inverted hierarchical spectrum in the case of a 
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Figure 4: Exclusion of an IH with mo = in the case of a true QD spectrum. We show contours 
for 1000 events (left) and 100 events (right) in the plane of the true Majorana phases assuming a true QD 
spectrum {m^^^ = 0.15 eV, sf^"^ = 0) fitted with IH and mo = 0, minimising with respect to all other 
parameters. 

true QD spectrum the branchings into tau leptons are crucial. If only electron and muon 
events are used in most regions of the parameter space an IH with mo = can fit data from 
a QD spectrum. For (ee), (e/i), (/i/i) branchings a degeneracy between IH and QD appears 
due to the freedom in adjusting S13, 6, 6*23 and the Majorana phases. This effect is also 
apparent from the black contour lines in Fig. [3] (lower-right panel). The significance of this 
degeneracy depends on details such as the errors imposed on s^g and S23, as well as on the 
systematical error cinorm- As discussed in Sec. 13.21 taking into account also decays into er 
and fiT is crucial to break this degeneracy, and in the full analysis used to calculate Figs. [2] 
and m the dependence on subtleties such as S13 and o"norm is small. 

3.4 Determination of Majorana phases 

Let us now investigate the tantalising possibility to determine the Majorana phases aij = 
ai—aj from the doubly charged Higgs decays. Since the decay is governed by a single diagram 
without any interference term the decays are CP conserving, and therefore no explicit CP 
violating effects can be observed. Nevertheless, the branchings depend (in a CP conserving 
way) on the phases, which eventually may allow to establish CP violating values for them. 
In general the measurement of Majorana phases is a very difficult task. Probably the only 
hope to access these phases will be neutrino-less double beta-decay in combination with an 
independent neutrino mass determination, where under very favourable circumstances [H] 
the phase ai2 might be measurable. 

We start by discussing some general properties of the branchings related to the Majorana 
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Figure 5: Determination of the Majorana phases for QD spectrum (toq = 0.15 eV) from 1000 doubly- 
charged Higgs pair events. We assume s^™" = and three example points for the true values of the Majorana 
phases given in each panel. The dashed lines in the middle panel correspond to the true values of the phases 
for which the degenerate solution according to Eq. (|4ip appears at a CP conserving value of 032- 



phases. Using Eqs. and 1^ one can write: 



BR^b oc \M, 



ab I 



(3^ 



From tliis expression it is evident that for a vanishing lightest neutrino mass, mo = 0, there 
is only one physical Majorana phase, 032 for NH and for IH. Next we note that since 
Ve3 oc Si3, it is clear that for S13 = all branchings involving electrons can only depend on 
C(i2^ Since the small effects of S13 cannot be explored efficiently, the determination of both 
phases simultaneously necessarily involves BR^^ and/or BR^^, see also Eqs. (!23l) to (1371) . 
Furthermore, from Eq. (138|) it can be seen that the branchings are invariant under 



a,; 



2-K — a 



(39) 



This symmetry is a consequence of the fact that there is no CP violation in the decays, 
and therefore the branchings have to be invariant under changing the signs of all phases 
simultaneously. 

In Fig. Elwe show that for a QD spectrum the observation of the decay of 1000 doubly- 
charged Higgs pairs allows to determine both Majorana phases. We assume some true 
values for the two phases and then perform a fit leaving all parameters free, where for S13 
we impose the constraint from present data. The actual accuracy to determine the phases 
depends on their true values, where we show three different examples in the three panels. 
For ai2 = 032 = vr (left panel) the allowed region is the largest, however the phases can 
be constrained to a unique region. In the other two cases the accuracy is better, but some 
ambiguities are left. The symmetry from Eq. 0391) is apparent in all panels, whereas in the 
case ai2 = 0:32 = vr it does not introduce an ambiguity. 



^For the same reason only ai2 can be tested in neutrino-less double beta-decay, where |Mee| is probed. 
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Figure 6: Determination of the Majorana phase for vanishing hghtest neutrino mass. We assume sf^° = 0. 
Left: 1, 2, 3(7 ranges for 032 as a function of its true value for NH assuming 1000 doubly-charged Higgs pair 
events. Right: 2, 3, Scr ranges for ai2 as a function of its true value for IH assuming 100 doubly-charged 
Higgs pair events. The dashed vertical lines indicate the region where CP violating values of Q12 can be 
established at 3a. 

The features of Fig. can be understood from Eqs. flHHl) to flH7|) . In addition to the 
symmetry Eq. (!39l) one finds that in the limit S13 = the phases 031 and 032 appear only in 
the particular combination 

2 

/ 2 2 \ /\ "^12 Sm Oii2 / . N 

[Si2 COS asi + cos 0:32) oc cos(a;32 — v^) with tan ip = 2 ; (40) 

^12 + "5l2 cosai2 

where we have used ayi and 0:32 as independent parameters. For constant a\2 there are two 
values of 0:32 which leave this combination invariant: for each 032 we expect a degenerate 
solution at 

«32 = 2v? - 032 • (41) 

For a\2 = 7r/2 one finds 2(f ^ 0.287r. In the case of 032 = vr shown in the right panel of 
Fig. [5] this degenerate solution appears at 0:32 — 1.287r, which cannot be resolved from the 
original one, and we are left with a two- fold ambiguity, due to Eq. flHIH) . In the middle panel, 
for ai2 = tt32 = 7r/2, the ambiguity (jHj) leads to a separated solution around 0:32 — I.TStt 
and, together with the symmetry from Eq. fl39p we end up with four degenerate solutions. 
However, in this case the individual regions are rather small, and the CP violating values of 
both phases can be estabhshed despite the presence of the four-fold ambiguity. 

Note that the symmetry (l39l) does not mix CP conserving and violating values of the 
phases, whereas this can happen for the degeneracy Eq. (jH]). The dashed curves in the 
middle panel of Fig. [5] correspond to the true values of the phases, for which 032 = or 
71. Hence, along these curves CP violating values for 032 cannot be established since the 
degeneracy is located at a CP conserving value. 
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Let us now discuss the potential to determine Majorana phases in case of hierarchical 
spectra. As mentioned above, in this case there is only one physical phase, a^2 for NH and 
ai2 for IH. In Fig. [6] we show the allowed interval for this phase which is obtained from 
the data as a function of its true value. In the fit the is minimised with respect to all 
other parameters. The left panel shows that for NH even with 1000 events at most a 2cr 
indication can be obtained, on whether 0:32 is closer to zero or tt. This can be understood 
from Eqs. (^Hj) to (ITTj) . which show that 032 appears at least suppressed by y/r. In contrast, 
as visible in the right panel, for IH a rather precise determination of 0:12 is possible already 
for 100 events, apart from the ambiguity ai2 — > 27r — ai2. For ai2 around 7r/2 or 37r/2 its 
CP violating value can be established, as marked by the vertical lines in Fig. El The good 
sensitivity is obvious from Eqs. ([2HD to which show a strong dependence of the leading 
terms in the branchings on 0:12. 

4 Summary and concluding remarks 

In this work we have adopted the assumptions that (i) neutrino masses are generated by the 
VEV of a Higgs triplet, (m) the doubly charged component of the triplet is light enough to 
be discovered at LHC, i.e., lighter than about 1 TeV, and {Hi) it decays with a significant 
fraction into like-sign lepton pairs. We have shown that under these assumptions LHC will 
provide very interesting information for neutrino physics. The reason is that the branching 
ratio of the doubly charged Higgs into like-sign leptons of flavour a and 6, BR(i7++ — > ^t^t)^ 
is proportional to the modulus of the corresponding element of the neutrino mass matrix, 
I Mat I ^. Hence the flavour composition of like-sign lepton events at LHC provides a direct 
test of the neutrino mass matrix. 

We have shown that the type of the neutrino mass spectrum (normal hierarchical, in- 
verted hierarchical, or quasi-degenerate) can be identified at the Scr level already with 100 
doubly charged Higgs pairs H H^^ decaying into four leptons. Typically such a number 
of events will be achieved for doubly charged scalar masses below 600 GeV and 100 fb~^ 
integrated luminosity, whereas for masses of 350 GeV of order 1000 events will be obtained. 
We have found that it is possible to decide whether the lightest neutrino mass is smaller 
or larger than roughly 0.01 eV, which marks the transition between hierarchical and quasi- 
degenerate spectra. If it is smaller the mass ordering (normal vs inverted) can be identified. 
A hierarchical spectrum with normal ordering has a distinct signature, namely a very small 
branching of the doubly charged Higgs decays into two electrons. Therefore, this mass pat- 
tern can easily be confirmed or ruled out at very high significance level. The other two 
possibilities for the neutrino mass spectrum, inverted hierarchical or quasi-degenerate, are 
somewhat more difficult to distinguish, but also in this case very good sensitivity is obtained, 
depending on the observed number of events. 

In this respect the inclusion of final states involving tau leptons is important, since if 
only electrons and muons are considered a degeneracy between IH and QD spectra appears. 
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In our analysis we have conservatively assumed that events where one of the four charged 
leptons is a tau can be reconstructed efficiently, thanks to the kinematic constraints and the 
information on the invariant mass of the event available from events without a tau. Certainly 
a more realistic study including detailed simulations and event reconstruction should confirm 
the assumptions which we have adopted here. 

The decay of the doubly charged Higgs in this framework does not show explicit CP 
violation, since the decay is dominated by a tree-level diagram without any interference 
term which could induce CP violation. Nevertheless, the CP conserving branching ratios 
strongly depend on the Majorana CP phases of the lepton mixing matrix. Therefore, the 
framework considered here opens the fascinating possibility to measure the Majorana phases 
in the neutrino mass matrix via CP even observables. Our results show that for an inverted 
hierarchical spectrum as well as for quasi-degenerate neutrinos this is indeed possible. In the 
first case, there is only one physical phase, ai2, which can be determined up to an ambiguity 
ai2 ^ 27r — ai2 already with 100 events. In the case of a quasi-degenerate spectrum both 
Majorana phases can be measured, where, depending on the actual values some ambiguities 
might occur. In many cases CP violating values of the phases can be established. 

Certainly the observation of a doubly charged scalar at LHC would be a great discovery 
of physics beyond the Standard Model. Of course this alone does by no means confirm the 
Higgs triplet mechanism for neutrino masses, since doubly charged particles decaying into 
leptons are predicted in many models. Therefore, in case such a particle is indeed found at 
LHC various consistency checks will have to be performed. It might turn out that the relation 
BR(if —>■ iaK) ^ l^aftP cannot be fulfilled for any neutrino mass matrix consistent with 
oscillation data. This would signal that a Higgs triplet cannot be the only source for neutrino 
masses. In this respect the information from decays into leptons of all fiavours (including 
taus) will be important. For example, also in the Zee-Babu model |1] for neutrino masses 
doubly charged scalars might be found at LHC. However, in this case branchings into tau 
leptons are suppressed by powers of (m^/m^)^ with respect to muons [27], whereas in the 
Higgs triplet model they are of similar size because of close to maximal 623 mixing. 

If LHC data on BR(iJ++ ^t^t) ^^^^ ^le consistent with a neutrino mass matrix from 
oscillation data, an analysis as pointed out in this work can be performed. Also in this case 
it will be of crucial importance to cross check the results with independent measurements, 
for example the determination of the neutrino mass ordering by oscillation experiments, or 
the measurement of the absolute neutrino mass in tritium beta-decay, neutrino-less double 
beta-decay or through cosmological observations. In particular, neutrino-less double beta- 
decay will provide a crucial test, since it gives an independent determination of the |Mee| 
element of the neutrino mass matrix, which — combined with information from oscillation 
experiments — will further constrain the allowed fiavour structure of the di-lepton events at 
LHC. The next generation of neutrino-less double beta-decay experiments is expected to 
probe the regime of the QD neutrino spectrum within a timescale comparable to the LHC 
measurement. Information from searches for lepton fiavour violating processes may be used 
as additional important consistency checks for the model. 
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In conclusion, a TeV scale Higgs triplet offers an appealing mechanism to provide mass 
to neutrinos, which can be directly tested at the LHC Such a scenario opens the possibility 
to measure the Majorana phases of the lepton mixing matrix, which in general is a very 
difficult — if not a hopeless task. 
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